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ABSTRACT 

This  paper  reports  Very  Large  Array  observations  at  325  and  1425  MHz  (A90  cm  and  A20  cm)  during  and  near  the 
periastron  passage  of  HD  80606b  on  HJD  2454424.86  (2007  November  20).  We  obtain  flux  density  limits  (3a)  of 
1 .7  mJy  and  48  /xJy  at  325  and  1425  MHz,  respectively,  equivalent  to  planetary  luminosity  limits  of  2.3  x  1024ergs  1 
and  2.7  x  1023  erg  s_1.  Unfortunately,  these  are  several  orders  of  magnitude  above  the  nominal  Jovian  value  (at 
40  MHz)  of  2  x  1 018  erg  s  1 .  The  motivation  for  these  observations  was  that  the  planetary  magneto  spheric  emission  is 
driven  by  a  stellar  wind-planetary  magnetosphere  interaction  so  that  the  planetary  luminosity  would  be  elevated  near 
periastron.  We  estimate  that,  near  periastron,  HD  80606b  might  be  as  much  as  3000  times  more  luminous  than  Jupiter. 
Recent  transit  observations  of  HD  80606b  provide  reasonably  stringent  constraints  on  the  planetary  mass  and  radius, 
and,  because  of  the  planet’s  highly  eccentric  orbit,  its  rotation  period  is  likely  to  be  “pseudo-synchronized”  to  its 
orbital  period,  allowing  a  robust  estimate  of  the  former.  Consequently,  we  are  able  to  make  relatively  robust  estimates 
of  the  emission  frequency  of  the  planetary  magneto  spheric  emission  and  find  it  to  be  around  60-90  MHz.  While  this 
is  too  low  for  our  reported  observations,  we  compare  HD  80606b  to  other  high-eccentricity  systems  and  assess  the 
detection  possibilities  for  both  near-term  and  more  distant  future  systems.  Of  the  known  high-eccentricity  planets, 
only  HD  80606b  is  likely  to  be  detectable,  as  the  others  (HD  20782B  and  HD  4113)  are  both  lower  mass  and  longer 
rotational  periods,  which  imply  weaker  magnetic  field  strengths.  We  find  that  both  the  forthcoming  “EVLA  low 
band”  system,  which  will  operate  as  low  as  65  MHz,  and  the  Low  Frequency  Array  may  be  able  to  improve  upon  our 
planetary  luminosity  limits  for  HD  80606b,  and  do  so  at  a  more  optimum  frequency.  If  the  low-frequency  component 
of  the  Square  Kilometre  Array  (SKA-lo)  and  a  future  lunar  radio  array  are  able  to  approach  their  thermal  noise 
limits,  they  should  be  able  to  detect  an  HD  80606b-like  planet,  unless  the  amount  by  which  the  planet’s  luminosity 
increases  is  substantially  less  than  the  factor  of  3000  that  we  estimate;  for  the  SKA-lo,  which  is  to  be  located  in  the 
southern  hemisphere,  future  planetary  surveys  will  have  to  find  southern  hemisphere  equivalents  of  HD  80606b. 

Key  words:  planetary  systems  -  planets  and  satellites:  magnetic  fields  -  planets  and  satellites:  individual 
(HD  80606b)  -  radio  continuum:  planetary  systems 


1.  INTRODUCTION 

HD  80606  is  a  G5  star,  at  a  distance  of  58.4  pc,  which  is 
distinguished  by  being  orbited  by  a  planet  with  one  of  the 
highest  known  orbital  eccentricities.  The  planet  is  a  giant  planet 
(3.94zb0.11M/)  with  a  1 1 1  day  orbital  period  on  an  extremely 
eccentric  orbit  ( e  =  0.93)  (Naef  et  al.  2001;  Pont  et  al.  2009). 

All  of  the  giant  planets  in  the  solar  system  and  the  Earth 
generate  radio  emission  as  a  result  of  an  interaction  between  the 
solar  wind  and  the  planetary  magnetosphere.  Specifically,  the 
solar  wind  impinging  on  the  planetary  magnetosphere  generates 
currents  within  the  magnetosphere  which  are  then  directed 
into  the  planet’s  magnetic  polar  regions  where  an  electron 
cyclotron  maser  instability  is  produced.  Motivated  by  the  radio 
emission  from  solar  system  planets,  there  have  been  a  number 
of  searches  for  magnetospherically  generated  radio  emission 
from  extrasolar  planets,  at  frequencies  ranging  from  25  MHz 
to  1400  MHz  (Yantis  et  al.  1977;  Winglee  et  al.  1986;  Bastian 
et  al.  2000;  Zarka  et  al.  2001;  Lazio  et  al.  2004;  Ryabov  et  al. 
2004;  George  &  Stevens  2007;  Lazio  &  Farrell  2007;  Smith 
et  al.  2009;  Lecavelier  Des  Etangs  et  al.  2009;  Lazio  et  al. 
2010).  Sensitivities  (la)  have  ranged  from  roughly  1000  mJy 
to  below  1  mJy,  with  generally  better  sensitivities  obtained  at 
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the  higher  frequencies.  Typically,  these  sensitivities  have  not 
been  sufficient  to  challenge  the  predicted  levels  of  extrasolar 
planetary  magneto  spheric  emission,  though  some  recent  multi¬ 
epoch  74  MHz  observations  of  t  Boo  (Lazio  &  Farrell  2007) 
do  indicate  that  its  planet’s  emission  is  not  consistent  with  the 
most  optimistic  predictions,  unless  the  emission  is  either  at  a 
frequency  well  below  74  MHz  or  is  much  more  strongly  beamed 
than  that  of  solar  system  planets. 

Most  recent  radio  observations  of  planets  have  focused  on 
planets  in  orbits  with  relatively  small  semi-major  axes  and  low 
eccentricities  (but  see  Lazio  et  al.  2010).  Part  of  the  motivation 
for  this  focus  is  the  stellar  wind  powering  of  the  electron 
cyclotron  maser  instability.  Close  to  the  host  star,  the  stellar 
wind  should  be  more  intense  with  a  concomitant  increase  in 
the  strength  of  the  planetary  radio  emission.  A  secondary  cause 
is  the  result  of  a  selection  effect.  Planets  close  to  their  host 
stars  produce  the  largest  radial  velocity  signature,  and  the  vast 
majority  of  known  extrasolar  planets  have  been  found  with  the 
radial  velocity  method. 

A  concern  with  this  focus  on  so-called  hot  Jupiters  is  that  they 
are  likely  to  be  tidally  locked  to  their  host  star,  with  rotation 
periods  of  order  a  few  days.  To  the  extent  that  rapid  rotation 
may  be  necessary  for  the  generation  of  a  strong  planetary 
magnetic  field,  these  “hot  Jupiters”  may  have  weak  fields  that 
do  not  generate  much  electron  cyclotron  maser  emission  or 


1929 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

DEC  2010 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 


4.  TITLE  AND  SUBTITLE 

Radio  Observations  Of  HD  80606  Near  Planetary  Periastron 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


6.  AUTHOR(S) 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 
5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

U.S.  Naval  Observatory  Flagstaff  Station, Flagstaff, AZ, 86001  report  number 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  Astronomical  Journal,  140:1929-1933, 2010  December 

14.  ABSTRACT 

This  paper  reports  Very  Large  Array  observations  at  325  and  1425  MHz  (&#955;90  cm  and  &#955;20  cm) 
during  and  near  the  periastron  passage  of  HD  80606b  on  HJD  2454424.86  (2007  November  20).  We  obtain 
flux  density  limits  (3&#963;)  of  1.7mjy  and  48&#956;Jy  at  325  and  1425  MHz,  respectively,  equivalent  to 
planetary  luminosity  limits  of  2.37024  erg  s&#8722;l  and  2.7  71023  erg  s&#8722;l.  Unfortunately,  these 
are  several  orders  of  magnitude  above  the  nominal  Jovian  value  (at  40  MHz)  of  27018  erg  s&#8722;l.  The 
motivation  for  these  observationswas  that  the  planetarymagnetospheric  emission  is  driven  by  a 
stellarwind?planetarymagnetosphere  interaction  so  that  the  planetary  luminositywould  be  elevated  near 
periastron. We  estimate  that,  near  periastron, HD80606b  might  be  as  much  as  3000  times  more  luminous 
than  Jupiter.  Recent  transit  observations  of  HD  80606b  provide  reasonably  stringent  constraints  on  the 
planetarymass  and  radius, and,  because  of  the  planet?s  highly  eccentric  orbit,  its  rotation  period  is  likely  to 
be  ?pseudo-synchronized?  to  its  orbital  period,  allowing  a  robust  estimate  of  the  former.  Consequently, we 
are  able  tomake  relatively  robust  estimates  of  the  emission  frequency  of  the  planetary  magnetospheric 
emission  and  find  it  to  be  around  60790  MHz.  While  this  is  too  low  for  our  reported  observations,  we 
compare  HD  80606b  to  other  high-eccentricity  systems  and  assess  the  detection  possibilities  for  both 
near-term  and  more  distant  future  systems. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Same  as 

5 

Report  (SAR) 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1930 


LAZIO  ET  AL. 


Vol.  140 


Table  1 

Observation  Log 


Frequency 

Wavelength 

Bandwidth 

Epoch 

Orbital  Phase 

325  MHz 

90  cm 

12.5  MHz 

HJD  2454423.89-2454424.22 
(2007  November  19  09:20-17:17  UT) 

0.69-0.73 

1425  MHz 

20  cm 

50  MHz 

HJD  2454424.74-2454425.07 
(2007  November  20  05:46-13:44  UT) 

0.94-0.12 

generate  it  at  a  frequency  well  below  that  accessible  from 
the  ground.  Absent  a  direct  radio  detection  of  magnetically 
generated  emission,  evidence  for  extrasolar  planetary  magnetic 
fields  is  ambiguous.  Consistent  with  the  scenario  of  weak  fields, 
France  et  al.  (2010)  were  unable  to  find  any  unambiguous  signal 
of  aurorally  generated  UV  emission  from  the  “hot  Jupiter” 
HD  209458b.  From  this  they  derive  a  limit  on  the  intrinsic 
planetary  magnetic  moment  of  0.01  AJj,  but  perhaps  as  large 
as  0.1  AJj,  if  the  interaction  with  the  stellar  magnetic  field  is 
important.  In  turn,  this  limit  implies  a  magnetic  field  strength 
of  only  10%  that  of  Jupiter.  Batygin  &  Stevenson  (2010)  argue 
that  Ohmic  dissipation  could  explain  the  inflated  atmospheres 
of  “hot  Jupiters,”  and  they  are  able  to  explain  the  atmospheres 
of  HD  209458b,  HD  189733b,  and  Tres-4b  assuming  a  field 
strength  of  10  G  at  the  cloud  tops,  as  compared  to  the  value  of 
4.2  G  for  Jupiter.  Strictly,  their  estimate  (~10  G)  is  an  upper 
limit  because  Ohmic  dissipation  need  not  be  the  only  effect 
operating  to  inflate  “hot  Jupiter”  atmospheres.  Lai  et  al.  (2010) 
show  that  the  asymmetric  transit  properties  of  WASP- 12b  might 
be  able  to  be  explained  by  a  magnetopause,  if  the  planetary 
magnetic  field  strength  is  of  order  1  G,  but  their  model  also 
has  a  fairly  weak  dependence  on  the  magnetic  field  strength 
(oc£1/3).  Conversely,  a  small  number  of  stars  with  “hot  Jupiters” 
show  evidence  for  star-planet  interactions  (e.g.,  Shkolnik  et  al. 
2008),  which  appear  to  have  a  magnetic  component.  In  contrast 
to  “hot  Jupiters,”  with  its  large  eccentricity  orbit,  HD  80606b 
may  be  immune  to  some  of  the  effects  that  might  act  to  damp 
the  magnetic  field  strength  of  “hot  Jupiters.” 

The  planet  orbiting  HD  80606  passed  through  periastron  on 
HJD  2454424.86  (2007  November  20  at  08:05  UT).  This  paper 
reports  observations  at  325  and  1425  MHz  (A.92  cm  and  A21  cm, 
respectively)  of  HD  80606b  close  to  its  periastron.  In  Section  2, 
we  describe  the  observations,  analyses,  and  present  the  images 
obtained;  in  Section  3,  we  discuss  our  upper  limits  in  the  context 
of  recent  multi- wavelength  observations  of  HD  80606b,  and  in 
Section  4  we  make  predictions  about  future  radio  observations. 

2.  OBSERVATIONS 

HD  80606  was  observed7  at  325  and  1425  MHz  with  the  Very 
Large  Array  (VLA),  in  its  B  configuration,  near  the  periastron 
passage  of  the  planet.  Table  1  summarizes  relevant  details  of  the 
observations,  with  specific  information  about  the  actual  times  of 
observations  and  the  planetary  phase.  In  addition  to  HD  80606, 
at  both  frequencies,  3C  286  was  observed  to  serve  as  a  flux 
density  and  bandpass  calibrator,  and  the  VLA  calibrator  J0929 
+  502  was  observed  as  a  visibility  phase  calibrator.  Although  the 
VLA  could  operate  at  74  MHz  (Kassim  et  al.  2007),  the  74  MHz 
dipoles  were  not  mounted  at  the  time  of  these  observations. 

Standard  data  reduction  techniques  were  used  within  the 
AIPS8  software  package  to  calibrate  the  spectral  bandpass, 
the  visibility  amplitudes,  and  visibility  phases.  The  data  were 

7  program  AL7 15. 

8  version  31 DEC08. 
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Figure  1.  Field  around  HD  80606  at  1425  MHz  on  HJD  2454425,  which 
includes  the  time  of  periastron  (Table  1).  The  cross  marks  the  location  of  the 
star.  The  beam  for  this  image  is  5749  x  4735,  the  noise  level  in  the  image  is 
16  /zJy  beam-1 ,  and  the  gray  scale  is  logarithmic  between  —80  /zJy  beam-1  and 
1.6  mJy  beam-1.  The  source  to  the  southeast  of  the  star  is  FIRST  J092239.6  + 
503529. 

acquired  in  a  spectral-line  mode  both  to  enable  the  identification 
and  excision  of  radio  frequency  interference  (RFI)  and  for  the 
purposes  of  wide-field  imaging.  These  issues  of  RFI  excision 
and  wide-field  imaging  are  particularly  relevant  for  the  325  MHz 
observations. 

Following  the  data  calibration,  several  iterations  of  imaging 
and  self-calibration  were  used  to  produce  images  of  the  field 
around  HD  80606;  at  both  frequencies,  polyhedral  imaging  was 
used  to  image  the  primary  beam  field  of  view  (Cornwell  &  Perley 
1992).  Figures  1  and  2  present  the  inner  portions  of  the  resulting 
images,  centered  on  HD  80606.  The  resulting  image  noise  levels 
are  16  /xJy  beam-1  at  1425  MHz  and  580  /xJy  beam-1  at 
325  MHz,  which  in  both  cases  are  within  approximately  50% 
of  the  expected  thermal  noise  levels. 

3.  DISCUSSION 

At  neither  frequency  do  we  detect  statistically  significant 
emission  from  the  direction  of  the  star.  Table  2  summarizes 
the  resulting  flux  density  limits  on  this  system  near  periastron 
passage  and  the  resulting  luminosity  limits.  For  reference,  the 
nominal  luminosity  of  Jupiter  is  2  x  1018  erg  s-1  at  a  frequency 
near  40  MHz. 

The  age  of  HD  80606  is  uncertain,  with  estimates  ranging 
from  as  little  as  0.3  Gyr  to  in  excess  of  14  Gyr  (Wright  et  al.  2004; 
Saffe  et  al.  2005;  Valenti  &  Fischer  2005;  Takeda  et  al.  2007; 
Holmberg  et  al.  2009),  but  with  many  estimates  favoring  an  age 


No.  6,  2010 


RADIO  OBSERVATIONS  OF  HD  80606  NEAR  PLANETARY  PERIASTRON 


1931 


09  22  55  50  45  40  35  30  25  20 


RIGHT  ASCENSION  (J2000) 

Figure  2.  Field  around  HD  80606  at  325  MHz  on  HJD  2454424,  just  before 
the  planet’s  periastron  passage  (Table  1).  The  cross  marks  the  location  of  the 
star.  The  beam  for  this  image  is  22'/ 1  x  20'/ 1  and  is  shown  in  the  lower  left 
corner;  the  noise  level  in  the  image  is  580  /idy  beam-1;  and  the  gray  scale  is 
logarithmic  between  —2.3  mJy  beam-1  and  58  mJy  beam-1 .  The  area  shown  is 
approximately  the  same  size  as  that  shown  in  Figure  1 . 


of  approximately  7  Gyr.  For  the  purposes  of  this  discussion, 
we  take  the  age  of  HD  80606  to  be  comparable  to  that  of  the 
Sun,  with  a  comparable  stellar  wind  strength.  Consequently,  we 
assume  that  a  Jovian  planet,  with  an  orbital  semi-major  axis 
comparable  to  that  of  Jupiter  (~5  AU),  would  have  a  radio 
luminosity  similar  to  that  of  Jupiter. 

It  is  clear  that  our  luminosity  limits  remain  several  orders 
of  magnitude  (105— 106)  above  the  luminosity  of  Jupiter.  The 
original  motivation  for  these  observations  was  the  expectation 
that  the  high  eccentricity  of  the  planet’s  orbit  could  produce  a 
dramatic  variation  in  the  strength  of  the  planet’s  radio  emission. 
Because  a  planet’s  magnetospheric  emission  is  controlled  by 
the  input  stellar  wind  power,  the  luminosity  of  the  planet  is 
a  strong  function  of  its  distance  from  the  star.  Farrell  et  al. 
(1999,  Appendix  A)  discuss  the  scaling  of  the  luminosity  with 
distance  and  find  that  an  approximate  scaling  law  is  dx ,  with 
x  ~  —1.6.  Griessmeier  et  al.  (2007a,  2007b)  consider  more 
complex  stellar  wind  models,  but  a  strong  dependence  on  the 
planet-star  distance  remains. 

Over  the  course  of  the  planet’s  orbit,  from  apastron  to  peri¬ 
astron,  we  therefore  expect  that  the  luminosity  of  HD  80606b 
could  vary  by  a  factor  of  200.  Even  at  apastron  (0.87  AU), 
HD  80606b  is  still  closer  to  HD  80606  than  Jupiter  is  to  the 
Sun.  If  we  consider  the  further  boost  that  could  be  obtained  by 
moving  a  Jovian-like  planet  from  a  distance  of  5  AU  to  0.87 
AU,  the  boost  would  be  a  factor  of  16.  Thus,  one  might  expect 
that,  near  periastron,  HD  80606b  could  be  as  much  as  a  factor 
of  3000  times  more  luminous  than  Jupiter.  However,  even  an  in¬ 
crease  by  a  factor  of  order  3000  over  Jupiter’s  luminosity  would 
be  insufficient  to  increase  HD  80606b  to  a  level  comparable  to 
our  detection  thresholds. 

A  simple  explanation  for  our  non-detections  is  therefore  that, 
even  with  the  added  boost  due  to  the  small  planet-star  distance 
during  periastron,  the  planet’s  luminosity  remains  too  low  to  be 
detected  with  current  equipment.  There  are  several  other  factors, 
however,  that  could  also  explain  our  non-detection. 


Table  2 

Luminosity  Limits  for  HD  80606b 

Wavelength 

Flux  Density  Limit 

Luminosity  Limit 

(3a) 

(erg  s-1) 

1425  MHz 

48  /zJy 

2.7  x  1023 

325  MHz 

1.7  mJy 

2.3  x  1024 

Notes.  In  converting  from  flux  density  to  luminosity,  an  emission  bandwidth  of 
order  the  observing  frequency  is  assumed. 


One  possibility  is  that  the  radio  emission  is  not  beamed  in  our 
direction.  In  the  case  of  Jupiter,  Zarka  et  al.  (2004)  determined 
that  it  emits  into  a  solid  angle  of  approximately  1.6  sr,  directed 
roughly  along  Jupiter’s  rotational  axis.  HD  80606b  is  one 
member  of  the  rapidly  growing  class  of  transiting  extrasolar 
planets,  and  recent  transit  observations  demonstrate  that  there 
is  a  significant  spin-orbit  misalignment  (~50°)  in  this  system 
(Fossey  et  al.  2009;  Garcia-Melendo  &  McCullough  2009; 
Moutou  et  al.  2009;  Pont  et  al.  2009;  Winn  et  al.  2009).  For  a 
transiting  planet  whose  rotational  axis  is  aligned  with  the  stellar 
rotation  axis  and  no  spin-orbit  misalignment,  the  probability 
of  intercepting  the  emission  beam  would  be  relatively  small. 
However,  the  existing  spin-orbit  misalignment  in  this  system 
increases  the  probability  that  emission  beam  intersects  the  line 
of  sight.  Of  course,  it  is  possible  that  the  emission  cone  of 
HD  80606b  is  much  smaller  than  that  of  Jupiter,  such  that  the 
probability  of  intercept  remains  quite  low. 

A  second  possibility  is  that  the  expected  upper  frequency  limit 
of  emission  is  (well)  below  our  observation  frequencies.  The 
intense  cyclotron  maser  emission  of  Jupiter  cuts  off  dramatically 
at  frequencies  above  approximately  40  MHz,  as  the  local 
plasma  frequency  in  the  polar  emitting  region  exceeds  the 
electron  cyclotron  frequency  and  the  radio  emission  cannot 
escape.  For  most  extrasolar  planets,  one  can  estimate  the  upper 
frequency  limit  by  scaling  their  magnetic  moments  from  that  of 
Jupiter.  Previous  estimates  for  the  magnetic  moment,  and  upper 
frequency  limit,  for  HD  80606b  have  varied  wildly.  Lazio  et  al. 
(2004)  estimated  an  upper  frequency  limit  of  around  1 80  MHz, 
while  Griessmeier  et  al.  (2007b)  estimated  0.8  MHz. 

The  combination  of  the  eccentric  orbit  and  the  planetary  tran¬ 
sits,  however,  now  provide  sufficient  information  to  estimate, 
with  relatively  low  uncertainty,  the  upper  emission  frequency. 
This  upper  emission  frequency  is  given  by  (Farrell  et  al.  1999; 
Lazio  et  al.  2004) 


vc  ^  24  MHz 


(1) 


for  a  given  planetary  mass  M,  rotation  and  radius  R.  In 
practice,  this  value  may  underestimate  the  actual  upper  emission 
frequency,  as  the  case  for  Jupiter  is  vc  ~  40  MHz  whereas 
Equation  (1)  suggests  24  MHz.  This  expression  also  assumes 
that  the  planetary  magnetic  dipole  is  not  offset  significantly  (as 
it  is  for  the  planets  Uranus  and  Neptune),  which  can  produce 
significant  asymmetric  emission. 

Of  the  various  determinations  of  the  mass  and  radius  of  the 
planet  (Fossey  et  al.  2009;  Garcia-Melendo  &  McCullough 
2009;  Moutou  et  al.  2009;  Pont  et  al.  2009;  Winn  et  al.  2009), 
we  adopt  those  that  Pont  et  al.  (2009)  obtained  from  a  Bayesian 
analysis  of  all  of  the  data  available  to  them.  They  found  a 
planetary  mass  of3.94±0.11  Mjandaradius  of0.98±0.03  Rj, 
implying  a  density  about  four  times  higher  than  that  of  Jupiter; 
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their  models  also  suggest  that  the  core  mass  of  HD  80606b  could 
be  several  times  that  of  Jupiter. 

Hut  (1981)  discusses  the  tidal  effects  in  eccentric  systems. 
The  large  eccentricity  for  the  orbit  of  HD  80606b  should  have 
driven  it  into  a  state  of  pseudo- synchronization  between  its 
rotation  and  orbital  periods.  Using  expressions  derived  by  Hut 
(1981),  we  estimate  a  (pseudo- synchronized)  rotational  period 
of  39.9  hr. 

Evaluating  Equation  (1),  we  find  that  the  upper  emission 
frequency  of  approximately  55  MHz,  and  potentially  as  high  as 
90  MHz,  if  one  accounts  for  the  50%  difference  between  the 
actual  value  for  Jupiter  and  that  predicted  by  Equation  (1).  It  is 
clear  that  this  upper  emission  frequency  is  well  below  that  of 
our  observations. 

This  estimate  for  the  emission  frequency  differs  from  that 
of  Lazio  et  al.  (2004)  by  only  a  factor  of  two  to  three,  but 
it  differs  from  that  obtained  by  Griessmeier  et  al.  (2007b)  by 
a  factor  of  order  102.  We  believe  that  the  differences  can  be 
attributed  to  assumptions  of  the  rotation  period.  Lazio  et  al. 
(2004)  assumed  a  rotation  period  equal  to  that  of  Jupiter,  which 
is  within  a  factor  of  four  of  the  rotation  period  we  adopt  here. 
However,  Griessmeier  et  al.  (2007b)  identified  the  planet  as 
“tidally  locked,”  which  would  mean  that  it  would  have  an 
extremely  slow  rotation  rate,  and  therefore  a  low  magnetic 
moment  and  low  emission  frequency.  One  of  the  key  factors  that 
may  have  led  Griessmeier  et  al.  (2007b)  to  identify  the  planet 
as  tidally  locked  is  that  they  appear  to  have  used  expressions 
for  estimating  the  synchronization  timescale  which  assume  a 
circular  orbit  (their  Appendix  B). 

4.  FUTURE  PROSPECTS 

We  consider  the  possibilities  of  future  observations  from  two 
perspectives.  First,  are  there  planets  other  than  HD  80606b  that 
should  be  targeted?  Second,  what  are  the  prospects  for  future 
instruments  to  detect  HD  80606b? 

While  HD  80606b  is  not  the  only  known  planet  with  a  high- 
eccentricity  orbit,  it  appears  the  most  promising  for  future  radio 
observations.  HD  20782b  is  a  1.9  Mj  (minimum  mass)  planet 
in  an  orbit  with  an  eccentricity  of  0.97  and  an  orbital  period  of 
592  days  (O’Toole  et  al.  2009),  and  HD  4113b  is  a  1.6  Mj 
(minimum  mass)  planet  in  an  orbit  with  an  eccentricity  of 
0.903  and  an  orbital  period  of  527  days  (Tamuz  et  al.  2008). 
The  pseudo-synchronization  rotational  period  is  coupled  to  the 
orbital  period,  so  these  longer  orbital  periods  imply  substantially 
larger  rotational  periods,  as  compared  to  HD  80606b.  We 
estimate  rotational  periods  of  63.8  hr  and  14  days  for  HD  20782b 
and  HD  4113b,  respectively.  Combined  with  their  lower  masses, 
Equation  (1)  suggests  that  their  emission  frequencies  are  likely 
to  be  below  20  MHz.  Unless  future  surveys  reveal  planets 
with  masses  and  orbital  eccentricities  comparable  to  those  of 
HD  80606b  (but  see  O’Toole  et  al.  2009),  it  appears  likely 
to  remain  the  most  promising  target  for  radio  observations  of 
(high-eccentricity)  planets. 

With  respect  to  observations  of  HD  80606b,  our  analysis 
suggests  that  future  observations  must  be  both  more  sensitive 
and  at  a  lower  frequency  than  those  that  we  report  here.  There 
are  several  observing  systems  that  are  either  under  construction 
or  in  the  design  phase  that  we  now  consider. 

The  Expanded  Very  Large  Array  (EVLA)  project  exposed 
certain  shortcomings  of  the  VLA’s  74  MHz  system  (Kassim 
et  al.  2007),  primarily  its  generation  of  harmonics  at  higher 
frequencies  that  fall  within  the  EVLA’s  operational  spectral 
window.  Consequently,  the  NRL  has  funded  the  National  Radio 


Astronomy  Observatory  to  develop  a  new  suite  of  low-frequency 
receivers.  With  the  previous  74  MHz  VLA  system,  automated 
data  processing  of  approximately  a  1  hr  observation  could  obtain 
rms  intensity  levels  of  order  100  mJy  beam-1  (Cohen  et  al. 
2007),  while  more  focused  efforts  on  longer  observations  could 
obtain  noise  levels  of  25  mJy  beam-1  (Cohen  et  al.  2004).  At 
a  minimum,  the  new  “EVLA-low  band”  system  will  have  a 
substantially  larger  bandwidth.  There  has  also  been  continuing 
algorithmic  improvements  in  the  areas  of  RFI  identification 
and  excision  and  ionospheric  calibration  (e.g.,  Athreya  2009; 
Interna  et  al.  2009;  Offringa  et  al.  2010).  As  an  ansatz,  we 
assume  that  a  combination  of  algorithmic  improvements  and 
the  larger  bandwidth  could  produce  an  approximate  factor  of  10 
improvement  in  sensitivity,  i.e.,  3  mJy  beam-1.  The  resulting 
(3cr)  luminosity  limit  would  then  be  1023  erg  s-1  (for  an  assumed 
10  MHz  intrinsic  emission  bandwidth  of  the  electron  cyclotron 
maser  emission).  This  value  is  comparable  to  our  more  stringent 
limit  (at  1425  MHz),  but  it  would  be  at  a  frequency  at  which 
emission  might  more  likely  be  detected. 

The  Low  Frequency  Array  (LOFAR  http://www.astron.nl/ 
radio-observatory/astronomers/lofar-astronomers)  is  currently 
under  construction  in  the  Netherlands  and  other  European 
countries.  When  complete,  it  will  consist  of  many  dipole  phased- 
array  stations,  operating  in  a  “high  band”  (HBA)  and  a  “low 
band”  (LBA).  LOFAR-LBA  band  will  cover  the  frequency 
range  30-80  MHz,  with  a  peak  sensitivity  around  60  MHz,  well 
matched  to  the  expected  emission  frequency  of  HD  80606b. 
Current  estimates  for  the  sensitivity  of  LOFAR-LBA  are  that 
it  will  produce  a  similar  or  slightly  better  thermal  noise  level 
(about  7  mJy  beam-1  in  a  1  hr  integration)  as  the  EVLA-low 
band  system.  We  emphasize  that  the  limiting  factor  in  imaging 
sensitivity  might  very  well  not  be  the  thermal  noise  level,  but 
other  factors  such  as  ionospheric  calibration.  Nonetheless,  it 
appears  that  the  LOFAR-LBA  band  system  would  produce 
similar  constraints. 

Following  the  EVLA-low  band,  LOFAR  and  other  pathfinder 
systems  are  likely  to  be  the  low-frequency  component  of 
the  Square  Kilometre  Array  (SKA-lo)  and  ultimately  a  low- 
frequency  array  on  the  Moon  (the  Lunar  Radio  Array).  The 
current  design  goal  for  SKA-lo  has  a  lower  operational  fre¬ 
quency  of  70  MHz  with  enough  collecting  area  that  the  notional 
thermal  noise  level  would  be  of  order  10  /zJy  beam-1  in  a  1  hr 
integration.  Unfortunately,  candidate  sites  for  the  SKA-lo  are  in 
the  southern  hemisphere  (Australia  and  South  Africa,  terrestrial 
latitude  ~  —  30°),  from  which  HD  80606b  cannot  be  observed. 
If  future  planetary  surveys  detect  other  HD  80606b-like  planets 
and  if  the  thermal  noise  level  can  be  approached  with  the  SKA- 
lo,  the  implied  luminosity  limit  would  then  be  3  x  1020  erg  s  ], 
well  within  the  range  we  estimate  might  be  plausible  based  on 
the  larger  stellar  wind  loading  of  HD  80606b’s  magnetosphere. 
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Inc.  The  LUNAR  consortium  is  funded  by  the  NASA  Lunar 
Science  Institute  (via  Cooperative  Agreement  NNA09DB30A) 
to  investigate  concepts  for  astrophysical  observatories  on  the 
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